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Abstract 


The  metnod  of  expansion  in  powers  of  a  thickness  , 

parameter  £  is  employed  to  calculate  supersonic,  in- 
viscid  flow  about  symmetrical  arrow-head  wings  lying  en¬ 
tirely  within  the  tip  Mach  cone.  In  this  method,  the  first- 
order  terms  in  €.  constitute  the  familiar,  linearized 

Prandtl-Glauert  approximation.  The  method  is  one  of  it era- 

* 

tion,  so  that  the  X  the  approximation  always  depends  on  the 
(X"^) th,  etc. 

Here  the  second  approximation,  i.e.,  the  terms 
in  £  and  ,  are  computed  and  plotted  for  a  family 
of  arrow-head  wings  having  various  leading-edge  angles  and 
thicknesses,  and  flying  at  various  Mach  numbers.  It  is 
necessary  to  use  Lighthill’s  method  to  determine  the  strength 
of  the  attached  conical  shock. 
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Introduction 


The  problem  to  be  considered  is  one  in  which  dissipative 
phenomena  appear  only  within  the  shock  wave.  Since  this  entropy 
change  is  of  third  order  in  thickness,  it  may  be  neglected  in  a 
second^order  theory.  Hence  the  entire  flow  field  may  be  considered 
isentropic  and  irrotational,  and  the  introduction  of  velocity 
potential  is  justified. 

The  exact  flow  over  a  conical  body  in  an  ideal,  uniform 
supersonic  flow  ?  .eld  is  conical.  By  expanding  the  solution  in 
powers  of  a  geometrical  parameter  £  ,  one  sees  that  each  term  of 
the  expansion  must  alee  be  conical. 

As  disturbances  may  not  propagate  upstream,  the  flow  over  an 
"arrow  head"  wing  is  identical  with  that  over  an  infinitely  extended 
body. 

That  the  solution  to  the  differential  equations  of  motion 
can  be  expanded  in  powers  of  thickness  parameter  can  be  justified  by- 
assuming  an  expansion  in  terms  of  diminishing  order  of  magnitude  and 
proceeding  to  satisfy  the  successive  boundary  conditions.  For  a  flat 
airfoil,  the  first-order  solution  is  formed  to  be  of  order  £  (thickness 
parameter).  This  is  dictated  by  the  boundary  condition  on  the  normal 
velocity  at  the  airfoil.  On  the  airfoil,  except  at  or  near  singular 
points  of  the  first  order  solution  (e.g.,  leading  edges),  the  same 
boundary  condition  then  requires  that  the  next  nontrivial  approximation 
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.be  of  order  <£  z  ,  and  so  on.  This  statement  may  be  checked 
against  the  boundary  condition  formulated  in  Eq*  (29).  Thus,  the 
power  expansion  in  £  is  the  only  expansion  which  will  satisfy  the 
boundary  conditions  on  the  body.  Recent  work  by  Lighthill  indicates 
that,  to  second  order,  the  power  series  expansion  in  £  properly 
describes  conditions  at  the  Mach  cone  as  well* 

Further,  it  is  expected  that  the  solution,  to  each  order  will 
be  analytic  at  some  plane  (e.g*,  a  plane  of  symmetry)  near  the  thin 
body;  thus,  the  value  of  each  approximation  on  the  body  can  be  ex¬ 
pressed  as  a  power  series  in  £  ,  the  leading  term  being  its  value  on 
the  plane  of  symmetry,  say.  As  will  be  shown  later,  one  may  then 
formulate  boundary  conditions  at  the  plane  of  symmetry  to  all  orders. 

Broderick  and  Lighthill  have  solved  corresponding  problems 
of  the  supersonic  flow  about  bodies  of  revolution  and,  in  order  to 
satisfy  boundary  conditions  at  the  surface,  have  found  it  necessary 
to  introduce  nonregular  terms  involving  the  logarithm  of  the  thickness 
parameter.  This  is  a  consequence  of  the  fact  that  the  axis  of  the 
body  (corresponding  to  the  plane  of  symmetry  in  the  present  case)  is 
a  singularity  of  their  solutions. 

In  the  case  of  supersonic  conical  flows  over  slender  bodies, 
an  attached  conical  shoek  wave  occurs  in  the  vicinity  of  the  leading 
Mach  cone.  In  the  linearized  theory  of  such  flows,  it  is  customaiy 

t 

to  assume  a  zero-strength  shock  located  at  Mach  cone.  Lighthill  has 
pointed  out  that  while  the  assumption  of  zero  shock  strength  is  valid 
for  the  first  (linearized)  approximation,  the  true  shock  strength  is 


WADC  TR  52-277 


2 


|s  second  order  in  disturbance  velocities.  He  presents  formulas  that 
provide  a  quantitative  second-order  description  of  the  flow  near  the 
shock  wave  in  terms  of  results  of  the  linearized  solution.  Lighthill * 
results  will  be  used  here  to  provide  boundary  conditions  at  the  Mach 
cone0 


The  pressure  coefficient  is  given  by 


c  = 


P-  Pc 


>G 


Pl  foQ  U 


=  itU)  + 


(l) 


Denoting  the  velocity  vector  by 

[is+ 7  *  M  [t  (2) 


where  the  unit  vector  ±  lies  in  the  streamwise  direction,  and 


f  =  /  M^-l  , 


(3) 


and  applying  the  isentropic  equations  of  motion,  we  have 


<r 


ci2>=  -U  (*%)  yYvi/,  (^rj 


(U) 

(5) 


Thus,  for  pressure  coefficient,  one  needs  in  addition  to  all  three 
first-order  cartesian  velocity  components,  the  second-order  correction 
to  the  streamnise  velocity  component  as  wello 


WADC  TR  52-277 


Notation 


? 

s 

GO 

5 

P 

fi 

C 

if 

M 

P 

M,  \Tj  W 

Z.k 

i 

(X- 

r 

\7X 

t,  r,  e 

l/1  /  P,  f>  A 
^  A  /C 


■  cartesian  space  coordinates 

■  conical  radial  coordinate 

=  Ts chaplygin  radial  cooordinate 

=*  space,  conical,  and  Tschaplygin  angular  ordinate 

■  complex  coordinates  in  the  Tschaplygin  plane 

■  pressure 

■  density 

*  pressure  coefficient 

*  free-stream  velocity 

*  Mach  Number 

■  a  constant  based  on  free-stream  Mach  Number 

»  cartesian  velocity  components  in  the  directions  ,  ,  and 

respectively 

■  angles  -  See  Fig.  2 

*  velocity  potential 

■  velocity  of  sound 

*  ratio  of  specific  heats 

*  the  Laplacian  operator 

»  geometrical  quantitites  in  the  Tschaplygin  plane  -  see  Fig.  3 
=  flow  functions 

■  constants 


Subscripts: 


oc.  signifies  evaluation  in  the  free  stream 

b  signifies  evaluation  on  the  bocfy 
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0  signifies  evaluation  on  the  plane  of  symmetry 

/>  signifies  particular  integral 

c  signifies  complemstary  function. 

The  subscript  notation  for  partial  differentiation  is  used 
•where  convenient. 

Superscripts: 

(nj  denotes  the  order  of  approximation. 

Primes  denote  ordinary  differentiation. 
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PART  I-THEOKT* 


A,  The  General  Case 

(l)  Differential  Equations  for  the  Velocity  Potential 

For  the  scheme  shown  in  Fig.  1,  the  isentropic  equations  of 
motion  are 


(cf-  9x2)  <?AA  *  (  <pp  +  (  a 

ft?  -  =  0 

^  1-  (  TJX  -  <px  l-  ^ (7) 


It  is  assumed  that 


f 


U*  +  £<pl°  +  t'f 


0r\j 


(8) 


Introducing  Eqs„  (7)  and  (8)  into  Eq.  (6)  and  collecting  terms  of 
order  £  and  £  z,  respectively,  one  finds 

/<?:  <n 


- 


XX 


CO 


=  o 


rxx  -  T-jy 
M*c  c> 


tAI  ('*%*)**'+  ?**♦</ 


(9) 


(10) 


Eq.  (9)  is  the  equation  of  the  linear-perturbation  theory. 

j\<i) 

The  solution  for  f  willbe  composed  of  a  particular  solution 
of  Eq.  (10)  and  a  complimentary  solution  satisfying  the  homogeneous 
counterpart  of  Eq0  (10),  which  solutions,  taken  together,  satisfy  the 
boundary  conditions  appropriate  to  the  particular  case  considered. 

?  7 
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(2)  Transformations  and  Particular  Integrals 

A  transformation  introduced  by  Tschaplygin  is  applicable 
to  conical  flows  and  changes  the  three-dimensional  hyperbolic  equations 
for  velocity  components  into  two-dimensional  Laplace's  equations*  The 
transformation  is  introduced  as  follows* 


J7'+z'’ 


60  =  > 


5  = 


r 


*+  j 


r 


(ii) 


Applying  this  transformation  to  Eqs*  (9)  and  (10),  it  is  found 
that 


V 

2 <LC>)  (  J. 

— 

= 

V x 

^‘°-o 

(12) 

,(Z)  2 
(7  aA.  ■ — 

fi  n- 

S 

■at;2 

r  +  s 

r/. 

■JVJ  4  at  7 

u)  _ 

/  I 

(13) 

( '~sV*  l 

^  /  1 

^  4s1- -sV 

C<r*3 

co 

+  S( !- 

UJ  Ass  - 

('Z- 

S  CO 

^  5  IAJ  J 

(Hi) 

u>=. 

( /*- 

sVaU 

CO  A.S 

] 

(15) 

4 

S  A+sly  LO  A  Jj 

(V 

-5  J  CUr?  co 

where 

4  H  A*  w"J"] 


(16) 


and 


/  + 


(17) 


The  following  transformation  may  now  be  introduced: 


J  =  se* 


'uJ 


S  = 


-  c'oO 
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(18) 


and  one  finds  ohat 


4  (  d  2 5 ) 


(19) 


■which  makes  it  possible  to  define  the  following  functions: 

ft  ( 5 J  v-  J,  (5)  =  OA  »">  , 

A  J  |  (20) 

f3  OJ  +  h  (5)  =  A  *r (n  ) 

and  to  write  a  particular  integral  for  •v.'1'*,  after  numerous  integra¬ 
tions  by  parts  in  Eq.  (13),  as  follows: 


Expressions  for  V~^  and  ^  ^similar  to  the  above  can  be  obtained 
from  Eqs<>  (lit.)  and  (15). 

It  is  clear  at  this  point  that  the  solution  given  in  Eq.  (21) 
will  be  infinite  at  the  tip  Mach  cone  -‘)  unless  the  derivative  in 
the  first  right-hand  term  vanishes  there.  This  is  an  example  of 
singular  behavior  that  concerned  Lighthill  in  the  investigation 
alreacty-  mentioned.  In  the  cases  considered  in  detail  here,  in  which 
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the  body  lies  entirely  within  the  tip  Ifach  cone,  this  derivative 
does  in  fact  vanish,  since  the  linearized  theory  leads  to  vanishing 
perturbation  velocities  at  this  cone.  Consequently,  the  second- 
order  solution  is  free  of  this  particular  difficulty,  and  one  needs 
Lighthill * s  technique  only  to  establish  the  upstream  boundary  condition. 

(3)  The  Boundary  Condition  on  the  Body 

One  prescribes  simply  that  the  flow  through  the  body  surface 
vanish  to  the  second  order  of  approximation.  It  would  be  convenient 
to  satisfy  this  boundary  condition  at  some  mean  plane,  as  is  custom¬ 
arily  done  in  the  linear  theory.  The  assumption  of  analyticity  of 
solutions  at  this  plane  enables  one  to  do  this.  The  details  of  this 
procedure  will  appear  later  in  an  illustrative  example. 

(li)  The  "Boundary  Condition"  at  the  Mach  Cone 

Since  it  is  intended  to  deal  herein  with  velocity  components 
rather  than  the  potential,  it  is  necessary  to  formulate  boundary 
conditions  on  these  quantities  at  the  Mach  cone.  As  has  alreacfy  been 
mentioned,  the  formulation  of  correct  boundary  conditions  to  be  satis¬ 
fied  by  the  perturbation  velocity  components  at  the  tip  Mach  cone  is 
based  upon  Lighthill^  detailed  analysis  of  the  conditions  near  the 
shock  wave.  It  may  be  of  value  to  review  briefly  the  salient  features 
of  his  method o 

The  correct  upstream  "boundary  conditions"  are  actually  the 
Mhock-wave  equations,  applied  at  the  shock-wave  location,  which  is 
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upstream  of  the  tip  Mach  cone  and  is  therefore  outside  the  region  of 
applicability  of  the  method  set  forth  above*  To  obviate  this  diffi¬ 
culty,  Lighthill  introduces  a  transformation  of  coordinates  from  y  ,  u> 
to  R  ,  w  ,  which  may  be  written  in  the  form 

>!=./?+  £  rc°  (  UJ)  y-  i  1  rtxJ(ui)  y-  •  •• 
where  the  successive  functions  Y n  Yu)  are  to  be  determinea  as  the 
calculation  progresses.  The  method  of  approximations  in  successive 
powers  of  <£  is  then  formulated  in  terms  of  the  new  independent 
variables  R  ,  u>  9  The  singularity  of  the  differential  equations  is 
thus  shifted  to  the  location  /?=/  ,  instead’ of  7  =/  *  The  choice 

of  the  functions  Ylni  ( to)  is  made,  however,  in  such  a  way  as  to  avoid 
divergence  of  the  method  at  ft  —  f  .It  is  then  found  that  the  shock 
wave  occurs  for  ft  -£  /  -  i.e.,  in  the  region  of  convergence. 

Lighthill  next  introduces  the  Rankine-Hugoniot  shock-wave 
equations  to  connect  the  conical  flow  to  the  flow  upstream  of  the 
shock.  This  leads  him  to  fictitious  boundary  values  assumed  by  the 
velocity  components  at  ^  =  /  •  In  particular,  consider  the  perturba¬ 
tion  velocity  potential 
y  -  U  y  f  (  ft, 

According  to  the  shock-wave  equations,  this  quantity  is  continuous 
across  the  shock  -  i.e.,  it  has  the  value  zero  just  behind  the  shock 
in  the  present  problem.  Lighthill  shows  that  the  corresponding  values 
to  be  applied  at  ft  - 1  are 

i‘°t  /,  w;  =  fUJ(/,co)  =  o 
fn*  (tj  =  0 
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Furthermore,  in  any  case  where  the  body  lies  entirely  within 
the  tip  Ikch  cone,  it  is  found  that  ( u>)  is  identically  zero. 

This  means  that,  at  any  point  in  the  flow  field,  "j  and  R  will 
differ  only  to  second  order  in  £  •  Thus,  away  from  the  vicinity 
of  the  shock  wave,  the  independent  variable  R,  in  the  functions 
may  be  directly  replaced  by  ^  without  introducing  any  error  greater 
than  order  . 

It  follows  from  the  fact  that  C<R)  is  zero  a  that 
the  differential  equations  satisfied  by  ~f'°  (R,  m)  and  lO) 

are  just  the  equations  .equivalent  to  Eqs.  (9)  and  (10)  and,  hence, 
to  Eqs.  (12)-(15)  (the  only  differences  being  those  due  to  the 
definition  of  the  variables).  Finally,  then,  to  the  order  of 
approximation  contemplated  in  this  paper,  it  will  be  correct  to  solve 
Eqs.  (l2)-(l5)  and  to  apply  at  f  =  /  the  boundary  conditions  (22). 

From  the  definition  of  /  '  J  (R,  co)  and  the  relation  (ll),  it 

follows  that  these  are  equivalent  to 
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Applying  the  Tschaplygin  transformation  to  the  irrotationality 
conditions 

~  ^  ^ 

and  evaluating  them  at  the  plane  cO  =  o  ,  rt  ,  one  obtains  the 

following  useful  relations  s 


(28) 
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B*  Arrow-Head  Inside  Mach  Cone 


(1)  Statement  of  Problem  and  Formulation  of  Boundary  Conditions 
Fig©  2  shows  the  configuration  of  the  bocfer  -  an”arrow-head" 
at  zero  angles  of  attack  and  yaw  with  respect  to  a  free  stream  of 
velocity  V  and  Mach  Number  'K*?  *  6  is  the  thickness  parameter. 

After  certain  geometrical  maneuvers,  the  boundary  condition 
at  the  body  for  /  >  °  may  be  written  as 

(Alt  =  ]>-  (  i<p,1 


7  JJ> 


(29) 


Assuming  analyticity  of  solutions  near  the  plane  of  symmetry,  one 
writes 


/  ?z  fv*#,2)]b  ~  t  in>]  +  e  ■+  t  V*  {x.pj 


(30) 


Substituting  expressions  (8)  and  (30)  into  Eq.  (29)  and  collecting 
terms  of  like  order  in  6 , 

l  A"  h  -  ^ 


tt”].-  -b+l  9 »:i.  -  ^a[^j 


(31) 


(32) 


Eqso  (31)  and  (32)  are  applicable  on  the  plane  of  symmetry  within  the 
leading  edges. 

Utilizing  the  Tschaplygin  transformation,  noting  that  5  =  / 
represents  the  Mach  cone,  and  defining  ^  as  the  S-  coordinate  of 
the  leading  edges,  Eqs.  (31)  and  (32)  may  be  written  as 
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(33) 


symmetty, 

(  ur'°]  a  =  I^J0  —<>  (  60  =  0'  L<s  <  I)  (35) 

xin)  is  even  in  both  "Z  and  £  n 

y is  even  in  z,  $  odd  in  £  r  (36) 

\AT!n)  is  odd  in  Z,  even  in  %  ' 

On  the  circle  s^/t 


The  differential  equations  applicable  in  the  Ts Chaplygin  plane 

are  given  by  Eqs.  (l2)-(l5) •  Use  will  be  found  for  the  irrotationality 

conditions  (21;)  through  (28)* 

Fig.  3  represents  the  Tschaplygin  plane,  in  which  the  problem 
will  be  solvedo  Stated  in  its  simplest  terms,  the  plan  of  attack 
win  be  as  follows: 

(a)  Find  ^  ^  and  Wu)3  which  satisfy  differential  Eq.  (12), 

and  the  irrotationality  conditions,  subject  to  the  boundary  conditions 
developed  above * 
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(b)  Find,  the  particular  integral  for  ^l'^9  using  the 
expression  (21).  To  this  will  be  added  the  complementary  function 
required  to  ensure  satisfaction  of  the  boundaiy  conditions,  thus 
completing  the  solution  for  M.(X\ 

(c)  Using  these  results,  apply  Eqs.  (It)  and  (f>)  to  find  the 
pressure  coefficient. 

( 2 )  Singularit ies 

Certain  difficulties  are  encountered  in  carrying  out  the 
above  procedure  because  of  the  presence  of  singularities  in  and 
and,  hence,  in  all  components  in  the  second  approximation,  at  the 
leading  edges.  This  is  dus  to  the  subsonic  nature  of  the  flow  at 
these  edges  and  the  consequent  stagnation  in  the  component  of  velocity 
in  the  plane  of  symmetry  and  normal  to  the  leading  edge.  A  perturba¬ 
tion  theory  always  fails  to  represent  stagnation  conditions  and  provides 
singularities  in  velocities  instead. 

One  would  expect  that  the  singularities  arising  in  this  problem 
at  the  leading  edtes  would  be  of  the  same  types  as  those  appearing 
in  the  case  of  an  infinite  yawed  wedge,  when  the  angle  of  yaw  is 
sufficient  to  provide  subsonic  normal  flow  at  the  leading  edge.  This 
normal  component  would  contain  the  only  singularities.  Further, 
the  subsonic  two-dimensional  flow  over  a  wedge  will  have  the  same  type 
of  singularity  as  that  of  the  infinite  yawed  wedge;  the  incompressible 
(/■7U  ==•  O  )  two-dimensional  flow  over  a  wedge  will  afford  the  same 
singularities  as  the  compressible  subsonic  case,  inasmuch  as  the  leading 
edge  in  the  latter  case  represents  a  stagnation  point  in  whose  neighbor¬ 
hood  the  local  Mach  Number  is  nearly  zero. 
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Examination  of  the  incompressible,  two-dimensional  flow  over 
a  wedge  shows  singjCLarities  at  the  leading  edge  of  the  logarithmic 
type  in  the  first  approximation  and  the  square  of  the  logarithm  and 
the  logarithm  in  the  second  approximation.  These  singularities  are 
all  integrable.  These,  then,  are  the  singularities  to  be  expected  in 
the  "arrow-head"  case0 

Admission  of  such  singularities  raises  the  question  of 
uniqueness.  It  is  known  that  a  potential  problem  has  a  unique 
solution  when  suitable  boundary  conditions  are  prescribed  on  a 
closed  contour,  provided  the  solution  is  to  be  regular  everywhere 
inside.  If  singularities  are  admitted,  an  ambiguity  may  be  expected  - 
i.e.,  certain  harmonic  functions  having  singularities  may  be  intro¬ 
duced  without  disturbing  the  boundary  conditions.  One  therefore 
inquires  what  harmonic  functions  having  singularities  of  the  type  jU/s-ij 
and/or  (  / 5  -  £  / )z~  on  the  axis  of  symmetry  may  be  added  to  the 

regular  solution  for  m‘  or  m'  ^satisfying  the'  appropriate  boundary 
conditions,  without  disturbing  these  boundary  conditions.  Thus,  in 
accordance  with  Eqs.  (24),  (36),  and  (37),  it  is  required  that  such 
a  singular  function  have  zero  angular  derivative  on  the  axis  of 
symmetry,  be  even  in  both  y  and  ,  and  vanish  on  the  unit  circle.  It 
may  be  shown  that  the  only  singular  function  meeting  these  requirements 
is  the  following: 

l<)  V  +  JU(r3)  \  ] 

(38) 
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Similar  reasoning  leads  to  the  following  ambiguous  term  for  xr<  ( or 

k '2  l  f  —  f  f  J  J  ( ^-  J  '^'"l  ^  ^  ^  y  (39) 

^  (n)  f*V 

The  constants  ^  and  /<2  have  superscripts  indicating  the  order 
of  approximation  involved# 

Since  w^&nd  vV^have  boundary  conditions  on  value,  they  can 
have  no  ambiguity  of  this  type.  Therefore,  knowing  iV,,;and  one 
may  find  A'//’  Ki^  and  /<J!J  .  Eqs.  (25)  and  (26)  are  con¬ 

venient  for  this  purpose. 

It  appears  that  the  lack  of  uniqueness  discussed  above  arises 
because  of  the  abandonment  of  the  velocity  potential  in  order  to  take 
advantage  of  the  simplicity  afforded  by  the  Ts Chaplygin  transformation. 

(3 )  First-Order  Solution 

(a)  Solution  for  <Ar'°-  The  differential  equation  is  given  in 
Eq.  (12),  and  the  boundary  conditions  are  given  in  Eqs.  (33) f  (35) ,  and 
(37)o  analogy  with  incompressible  fluid  theory,  the  solution  may  be 
written  down  immediately,  using  the  stream  function  for  the  incomjr  essible 
source  1 

vfU>  =  ~  (  V/n)  (  &•  •*■&>-  i-  9  3  +  ~  7T) 

7  m 

(b)  Solution  for  it10-  The  differential  equation  is  given  in 
Eq0  (12)  above.  The  boundary  conditions  qre  those  given  in  Eqs.  (25), 
(36),  and  (37),  and  from  Eq.  (2J4.)  s 

=  0  (  s  <J)  (ia) 
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It  is  plain  that  the  solution  analytic  inside  the  unit  circle  is 
identically  zero,  leaving  only  the  constant  Kl,'\o  be  found.  Applying 
Eq.  (25)  to  Eqs.  (38)  and  (Uo),  and,  for  convenience,  evaluating 
terms  near  the  leading  edge,  one  may  iramedL  ately  -write 

^  pn  ^  (irj  -  (try)]  (12) 

(c)  Solution  for  Qr  a  procedure  similar  to  the  one  employed 

to  find  it  is  found  that 

(to) 

(d)  Information  Obtained  from  the  Above  Solutions.  -  Expanding  (liO), 
one  finds  that 

and  from  Eq.  (Ii2)  is  obtained 


Eq.  (It5)  may  be  substituted  into  the  results  of  Lighthill  for  the 
strength  and  location  of  the  shock  -nave. 

For  use  in  finding  U  ^  ,  one  puts  /u°9  and  nT^into 

complex  form  to  conform  to  the  expressions  (20^  above. 
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(1;)  Second-Order  Solution 

In  philosophy,  the  solution  for  /J^is  obtained  in  the  same 
way  as  was  jUu'\ in  the  preceding  paragraph.  In  detail,  the  procedure 
is  as  follows*  One  may  write 

^  -h  ) 

«rn>  =  I  <*> 

IZ)  .  t "M 

where  is  obtained  from  Eq.  (21) ,  wp  is  obtained  from  a  similar 
expression  and  iind  vfl  satisfy  Laplace’s  equation  in  the  Ts chaplygin 
plane. 

(a)  Find  M-j? ;  note  that  the  functions  obtained  from  the 
integrations  indicated  in  Eq.  (21)  may  be  -complex.  i/f  .  may  be  found 
in  real  form  fcy  adding  to  these  complex  functions  any  convenient 
functions  of S  or J  alone. 

(b)  Find  ^t^such  that  ^  ^  ^  satisfies  the  boundary 

conditions  (21*),  (36),  and  (37)  and  has  no  singularities  of  order 
higher  than  (  /s-  ?!  )  z  .  Expression  (38)  is  to  be  included 

in  M  l'>. 

(c)  Find  J  . 

(d)  Find  vf \  ^  such  that  ^  ■+  W~t  satisfies  the  boundary 
conditions  (3k)  f  (35) ,  and  (36)  and  has  no  singularities  of  order 
higher  than  (  JL^  t l/) 1  .  Note  that  no  boundary  condition  on  W 

at  the  unit  circle  is  applied. 

(•)  Apply  Eq.  (26)  to  find  /C/  i  This  step  is  simplified 
by  considering  only  those  terms  that  are  singular  at  the  leading 
edge.  The  equation  will  contain  terms  singular  at  S  =  Jl  ,  and  one 
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may  require  the  equation  to  be  satisfied  upon  arbitrarily  close  ap¬ 
proach  to  the  point  JT  =  X  •  -^11  terms  nonsingular  at  this  point 
may  then  be  dropped,  and  all  nonsingular  factors  in  singular  terms 
may  be  evaluated  at  J  =  i  «  ^he  term  involving  the  unknown  K]ZJ 
has  a  derivative  with  a  simple-pole  type  of  singularity.  Thus,  one 
needs  only  to  formulate  Eq.  (26)  for  terms  having  simple  poles  in 
derivative.  Finding  /c/^in  this  way  makes  it  possible  to  ignore 
ary  boundary  condition  on  ur^at  the  unit  circle,  since  its 
application  would  result  in  a  regular  function  that  would  make  no 
contribution  to  the  determination  of  ^J. 

This  procedure  has  been  carried  out  and  formulas  have  been 
obtained  from  which  computations  can  be  made  (Part  Il)« 
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PART  II  COMPUTATION  AND  RESULTS 


A0  Computation  and  Results 
In  the  cases  computed,  the  following  values  are  used* 

u-t,  y^/,4- 

ST,  cJ. 

L  ~  a.  / i  0.3  ,  o.  78  . 

<3  L  -  6  a,  >o° . 


Introducing  the  notations 

/  -  rSLSSbtS 

V  IT  (5 

srf,j  _  us  n-^ 

'  -2  ~  “  if  ~7— 


n  /. 

'  fSn  i-z1 


/7  <5, 

/  2 

^3 

/  r1-  / 

f! 

-  / 

-/ 

2  *  / 

-  / 

~/ 

•v-  / 

J  4-/ 

-/ 

*•  / 

-  / 

2  V  I-  l1-  *  2n  C 


A*  =  tZEr,2.  Al-ieff,*  f  '  -f  .<>  .  f  -  ) 

As-^^ra.  t>£  *<>4,  »  f«y^O  .  ^  -  ,-TTi  Z 

^  ^  ^  tL  -  Ar) 


^  S’  >tJlK  / 

y  l  j 

=  ~c~t^SL  £n  Zj 


2iv  3  z 

f,+  l^Z.En  <£*? 


\  =  -  j-  tt-tvi  eSf  i+hi  -ia  <£  +  ) 

•C  /  12  -/W  V 


l£Sij 
r~>  —  ->  -I— / 1  4  r  1 


P>  =  -  i  ±±X  ?  Fx  l 


p+  2  i  Z  ( -+ zs  +  z? 
pjr  =  ~  '-r-l  =(>-*')'£  E* 

*.'[<+  *sO~^o] 

D*=~  Ds=jr%E~(>+sj"-f- 

o  ~  ~2*i  ES  (  >t  *,/-  Cj 
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and  farther  putting  the  form 

Mr'>  =  +  9*3 

the  procedure  of  computation  can  be  summarized  as  follows* 

a)  Evaluate  )  0  £°r  a  series  of  values  of  s  according 
to  the  formula  (p.  llil;,  Ref.  2): 

■4-f#),  =n>  ^.[LU^  +  Utr^-Unj- //. <h -  L<i>) 

+  V’«-7 ZTxlJ*  (Slj  -l~(Sr+)]£  ny  +  J*'rV] 

+  s[/»rx-£,rj]J  t  '(LiUnJ-i^ftr+j]] 

+  nb  j Iji  /►»  (-Pnj][ *  l-ej  M'ijJl—  ^  tj] J 

~f  5  Xy%  ( fl_)  -^v,  4  C n^)  -t-  *  <f  1 5  "*■  a  4"*  f  ^ 

+  SJL,  f0)Jt  [n/o  t~(4r,)  +  o9  +  *ff  ^rrjJ  +  &*'*'*)] 

+  o,,  jC 4, a  - l^Unj] l+  {&*  0  -4,  <^/vo7  V 

»/:i  /  a,  A  -  JUftr+y  J/"  4,0  - 

+  C4,rx  -  UGl  -L.f/r+1]  +  *,+  [&««,- Q  J  +  o,r  [6*  ft,  -<LRx 

+  nu,  /4,  r^i-  ly,o]j  jLUn)  +t*{*r>t-)  +  {&***_+  £4,0  -  2(^Q,  +  Ai  /?0j 

*  y'  /«?  -  ^f^fjtr.J  *~*~(  lr«)  *  4  A  *4, 0  7  j 

x  *j,7/,/hj7/'4/^J^4.a#./^o  -j£i^  -2[A.Q,+ 

•*  c-l.  [  4,  Q ,  +  4,  /?j.  J  i"  a  ty)a 
■h  £  k(,X)[ 4,  A  +•  4.  O  -  t~,Ur,J  - 

Vj  kjt-r-*. 


-  4, 2 

„  '+1^/  F  x  J<L 

n3  S-A4  ”,s  =  ,171  r-  *jr  /T7 »• 

=  £tl 

„  _ 

*  ~  f-g.  * 

o?  =  As~  nlt  _  _  a  £2/  '*tz 

n'°  =  _  - 

^  ^  ~  ^  n  '  b 

V=  -£** 

O,,  =.  A  7 

>1  r  »•  ^ 

^  ,-4 
**  =  4/ 

n ,  2.  As 

j£  /  —  2-  £  \ 

^/3  S  -  /4^  t  7~u  «  "  4/?7i/ 

"  V  =' 

>0  _  J.  -— -  /  c  '•  J  4  _  -^x  \ 

—  z  t~  TZJi-  L8*  / 

UADC  TR  52-277 


■  23 


(<P).  =^I.T7T?.7aV -4  *'"] 

I 

c,  -  -  tiej] 

+  ^ f*+ *'j  ~  ] 

-  A,  -j4ji  t  Aft  +  Ei'yZj*.  ^  c* I) 

^i.  £J  2  *~(/-  ^  V) 1  £  ^  f  '*•* )  ~  ( 1 

1-  ~  flji/C  §-  pei.j^o+*)-(4'+/f£^e~0-'ej] 

k,  -£  =  A3-£iV  -r  Air  (l^f  ■>-£)-  {]  +  At  -/U-— ^  ■+  Z-^yA  n~l  \) 
tA-g  fA,  (>-/'’}  -A*,  fi.4j ]  —  Ay  [ A<n( -A-,n—t'L j  ] 

+0^'  -'¥^'£rj] 

*  'ir-[z§^^( Ei)  '-*-  .  -  /ttl  TZr  ^'J 

2/^  (  z 

f J  ^  ^  f  /  “  -/  *^)  —  ^Vv  /*  'j 

~£i-'~ [£,(/+- /y  - £~r/~  j'j )]  } 

+  -^-(  £x[ -A _ £_-_^  )  _  £v_^  *-  _ '-// 

i-<L'-'c>li+X'-  ,-jx  xjj  c*Lu{>-  /-r  W  7/ 


■Z±£  ^  /£-  E  i  \ /_n  t3 

d-iy'  C  ‘  2)  * 


J-L  +  An  *  n  /X- 

A-iy*  +  +  * 


+  /-/*• 
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b)  Tabulate  (&)0  Tor  series  of  values  of  oo  from  the 
formula  (  p.  1U8,  Ref*  2)x 

~  4( fx)6  =  v  *».  /  &2l  +  frsJfa] 

+  4,  (*■)£  tk^rL  +  ^  o]  f  JL,  tQ)g[ *9 

+  */0  [ £,$]*+  *„  r^*  % J^+ 

-t  *»j  ^  -4o  7^  -  -^'5  J  *  *'*  {"-0L-Q]  +  f*,rt^ri  +  l'rj] 

t  Ay6f^-%7‘'-*  ^(Q,)Uk,*i  +( tbr-  *"i)^G] 

4  **,&-  £,7J4h-0.  +  c &,$*-  hi, 3 

i-  c±t  l^R,  r-^R^  ]  4  bii 

.  ■l  .  jc  n~f)T-  1 


m,  s  -  r*  +A>i 

A u  =  af  j Lv-  <4- / 

*j  =• 

<k(,  H  A-l.  /  %  1 

hn  =  -  +  AjJx) 

kg  =  -  (A2  +A  4) 


*/+-  -  71  f  FJ(r-Wx  f>  /-4  u  ) 

u  ,  —  _  '  >+ 1  u 

rrj\ 

-  _  />  /A'v  .  '_£iA  7 

At?  =  LtA? 

J-  TT  t~4  ’>- 

=  -ZE* 


*?«  s  /4i —  A  fc 


^/o  s  -  (Afii-Ai  +  E/fL 

h  ~  r r=  i_ii  /  £/  j- .7 

''  “  Z  /,«^vjV-  *  45V,  /-/wJ 
*'*■  2  =  *«,. 

Then  since  f is  even  about  both  vertical  and  horizontal  axes, 
^oe  may  be  expressed  in  a  Fourier  series 


-  2  4  •  4-  <E-  Ql^  C&3.  3  cO 


where 


Gtr^  <§l  *t  6<J  o(  CO  — 


=  ^ 


C«  doJ 
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is  evaluated  by  graphical  integration  using  a  planimeter. 

Thus  the  regular  harmonic  function  associated  with  boundary  value  (l) 
is 

=  2  4  b  +■  2[.  4  -I  S  *  ^  C-ir2  2  m  lO 
/ 

o< 

where  ( #  =  {  4  e  -r 

(A),  is  then  evaluated  at  the  series  of  values  -s  picked  above# 

«)  A  Fourier  analysis  of  is  made  for  Lighthill’s  result 
on  Mach  cone  (Appendix),  which  gives  on  ^  =  o  , 

C&0  =  5M0  y  £ 

where  /inis  obtained  from  formulas  furnished  in  Appendix. 

d)  Results  of  a),  b),  and  c)  above  are  added  to  give  ( m.U' )0 

e)  ^  and  oMare  computed  from  (I|2),  (JU3)  • 

f)  cf>  >  and  cp  are  computed  from  (1;),  (5) ® 

In  performing  the  calculations,  resort  has  been  made  to 
various  geometric  relations  furnished  by  Appendix  A,  Ref*  2. 

The  computation  procedure  is  illustrated  by  tables  for  the 
case  of  X  =  o-  3  .  The  computed  results  for  (A")  9 

Cf/'  *9  C for  various  Mach  numbers  and  values  of  A  are  plotted 
in  Curve  1  to  Curve  20.  Note  that  different  scales  have  been 
adopted  in  plotting  the  curves.  For  ,  plots  of  ^  and  Cp 

at  5  =  o  against  L  are  shown  in  Rig.  21.  For  X  -  o  3  , 

plots  of  k-  and  cya  at  *=o  against  are  given  in  Fig.  22.  In 

both  cases  2  t  =  6  °  is  used. 
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B*  Concluding  Remarks 


The  theory  and  ‘  computation  formulas  adopted  in  the  present 
report  are  mainly  taken  from  Ref,  2,  by  Ur,  F.  K.  Moore.  Those  not 
to  be  found  there  are  given  in  the  Appendix  to  the  present  report. 

The  boundary  condition  on  Mach  cone  given  in  Ref,  2  is  incorrect. 
Prof,  W,  R.  Sears  pointed  out  that  Lighthill 1  s  result  in  Ref,  3  should 
be  introduced.  Ref,  1  is  the  corrected  version. 

The  computations  have  largely  been  carried  out  by  Mrs,  Anne 

Kane, 

The  effect  of  <£  on  the  difference  between  the  first  and 
second  approximation,  as  expected,  grows  with  £ 

The  reason  that  the  discrepancy  between  first  and  second 
approximation  changes  radically  with  decreasing  planform  angle  A 
is  that  the  parameter  £  is  not,  in  fact,  an  ideal  thickness  parameter, 

A  little  r  eflection  shows  that  <£  '  =  £  A  might  be  a  better 
choice.  Indeed,  the  boundary  value  on  the  singular  line  inside  the 
unit  circle  of  the  conical  plane  is  given  by  £ 'rather  than  2  »  Thus 
even  fihr  such  small  values  of  L  as  might  seem  evidently  to  meet  the 
requirement  of  linearization,  the  value  of  £'  can  still  cause  break¬ 
down,  provided  A  be  sufficiently  small.  It  is  suggested  therefore 
that  the  present  second  approximation  might  be  improved  by  expanding, 
instead,  in  powers  of  £, 

It  might  also  be  of  interest  to  note  that,  so  far  as  computdion 
results  indicate,  the,  present  theory  predicts  a  particular  value  of 
■l  (  ,==  c.  3 )  ,  for  which  the  first-order  solution  gives  good  approxima¬ 

tion  irrespective  of  Mach  number,  A  further  investigation  to  s  ee 

whether  this  is  a  coincidence  is  also  suggested, 
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APPENDIX 


In  order  to  utilize  Lighthlll ♦ s  result,  it  should  be  noted 
that  in  his  discussion,  physical  polar  coordinates  ( J",  u>)  were  used, 
while  in  the  present  report,  conical  coordinates  ( .s,w)are  used.  The 
problem  can  be  stated  in  following  manner.  In  the  differential  equation 


^  2  /U.  J  =  /V  ( r,  lo) 


(1) 


a  solution  of  the  form 


^  =  t  +  9?  fs 


is  sought,  where 


V  2  P j  ~  o 


(2) 


(3) 


and  satisfies  the  boundary  conditions 

~  -  Ct  (  lO) 


($1 


( oo~  q) 


-  o 


(U) 

(5) 


Q  (-<*>)  is  obtained  from  Lighthill’s  result  on  the  study  of  conical 
shock.  Thus,  since  p  is  analytic  inside  the  unit  circle  *-/  , 
it  takes  on  the  form 


c< 


i>2  =  2  3" 


C-tr^  -?  oO 


(6) 


Now  using  Lighthill’s  notation,  Refo  3,  we  have 

(*■)  _  -i  r  f  J.  „  J-  \ 


~  ~LT  (  ~F  —  n  f^) 


u'  '  (si  =/J  -  r  IS  n.  /3.  -  -  <R-I) 
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(8) 


s°  ( £  ) 


%  -t-'  M  H  N  Cr^cO) 

2.  /?->-  A  ~  U 


1  over  r  =  -s  ■=■  J 

Bit  as  can  easily  be  shown  (Ref*  3) 

A  (  co)  =  (  =r^  J 

^  J  t-n.  s  ^  =/ 


and  since 


"7^n-  pi X  10  +  uj)  —  ( \J  Co-*  u j  1-  vA/UA,  coj 

r  **L.\  _  *± JE^ri  .  .  ^ 

j  - - J  —  t/  £  CaJ  - -  1 

^//Ta'/l*-/  TT  <7,  / 


/  -4*1  \ 

^7? 

^v/7-a  7l-/ 

7T 

,  0(V 

{  4 / i-fC  « / 

.  7> 

-/w« )  =  +  +£2^^) 

q  -J  l-n-  y^-i  pV  '■<V/-'v  <J]  t~  n.  J 


4JT  /  ■+  _ 

n  C<rxi'tO 


1,e*  C7feoJ  =r  -Y-'^Va))  = 


&+/)/£,  _ 

-"Y*  owy-*<w  <*>?«] 


~  <>_ /4in  °cr2  *■ 
6 

and  finally  we  arrive  at 


u 

where  4  ^  <'''",>  <"  ~AH 

/s  _  iL  ry.  .,  /_^»  y  v _ 

2-  r  )  i  pJ  (/-(*)  1  (V-/ /; 

/4yu  =  zi_rd>,u*y.  Itbt3-iy 
W"  0+ ){yj  (/-*'}*  ft- Jt*) 

Ah  =, 

^  7y.  </-  p'-p v  <7- 1  vj 

/4«p  =  JJz-fg-,  ■)  .Htf  ^  rr-3ly+) 
•7T'-(°*J  'pJ  Cf- 1  y  »■  t>~  -t  V 
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(18) 


a  (A.  /*-  w AL\<  1 
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( /-  A  Off*  U  J  i 
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(i  ”?>)(- >h^)<n 


_  _  2  *■  X'  (  ?  rn~hls-  /)  d?  _  _  2  i'  '  X'  iD 

—  ~a~i  (  ^jl  o  -+  j  ) 


ss  =/ 


JT-.S  =  - 


jr  = 

2Vy-  Ct  i- 


/  -  \J_i~ 


S  -hS  - 


/  ^  v  /-  O 
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